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Geo-neutrino studies are based on theoretical estimates of geo-neutrino spectra. We propose a 
method for a direct measurement of the energy distribution of antineutrinos from decays of long- 
lived radioactive isotopes. We present preliminary results for the geo-neutrinos from 214 Bi decay, 
a process which accounts for about one half of the total geo-neutrino signal. The feeding proba- 
bility of the lowest state of 214 Bi — the most important for geo-neutrino signal — is found to be 
po = 0.177 ± 0.004 (stat) qo? ( svs )> under the hypothesis of Universal Neutrino Spectrum Shape 
(UNSS). This value is consistent with the (indirect) estimate of the Table of Isotopes (Tol). We 
show that achievable larger statistics and reduction of systematics should allow to test possible 
distortions of the neutrino spectrum from that predicted using the UNSS hypothesis. Implications 
on the geo-neutrino signal are discussed. 

PACS numbers: 13.15.+g, 14.60.Pq, 23.40.Bw, 92.20.Td, 87.56.bg 
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I. INTRODUCTION 

Geo-neutrinos, the antineutrinos from the progenies of 
U, Th and 40 K decays in the Earth, bring to the surface 
information from the whole planet, concerning its content 
of radioactive elements. Their detection can shed light 
on the sources of the terrestrial heat flow, on the present 
composition, and on the origins of the Earth. 

Although geo-neutrinos were conceived very long ago, 
only recently they have been considered seriously as a 
new probe of our planet interior, as a consequence of two 
fundamental advances that occurred in the last few years: 
the development of large extremely-low-background neu- 
trino detectors and the progress on understanding neu- 
trino propagation. From the theoretical point of view, 
the links between the geo-neutrino signal and the to- 
tal amount of natural radioactivity in the Earth have 
been analyzed by several groups. Various "reference 
models" [3, 0, 0| for geo-neutrino production have been 
presented in the literature; some of these models have 
been refined with geological and geochemical studies of 
the regions surrounding the detectors Q. KamLAND 
and Borexino [7, 8] are collecting geo-neutrino data, 
while several planned experiments (e.g., SNO+, LENA, 
HANOHANO, EARTH, . . . ) have geo-neutrino measure- 
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ments among their primary goals. A recent review is 
presented in [9| ■ 

This activity has to be complemented with some deep- 
ening of the nuclear physics that is at the basis of geo- 
neutrino detection and which is crucial for interpreting 
future geo-neutrino data. The aim of this paper is to 
discuss the uncertainties of some nuclear physics param- 
eters that enter in the interpretation of the geo-neutrino 
signal and to provide a framework for an experimental 
determination of these parameters. 

In all experiments that use hydrocarbons as detection 
media, either running or in preparation, the reaction for 
geo-neutrino detection is the inverse beta decay on free 
protons 

u e +p^e + + n- 1.806 MeV . (1) 

The signal is estimated from the cross section ct(_E„) 
of Eq. |T]) and from the decay spectra f(E„) of geo- 
neutrinos produced in each beta decay along the decay 
chains, the relevant quantity being the "specific signals" 
defined as: 

Si = f dE vt T(E v )fi(E v ) , (2) 

J Eq 

where Eq = 1.806 MeV is the threshold energy for reac- 
tion (JTJ) , E max is the maximal geo-neutrino energy and 
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the spectrum f(E v ) is normalized to one geo-neutrino 1 

/ dE v f{E v ) = \ . (3) 
Jo 

It is important to observe that the specific signal is af- 
fected by unknown uncertainties. In fact, whereas o~(EjA 
is affected by uncertainties of less than percent [Tol. 
it is difficult to assess the accuracy of f(E v ), which is 
determined from rather indirect measurements and ques- 
tionable theoretical assumptions. 

Our goal is to provide a framework for a direct mea- 
surement of f(E v ), so that the accuracy of the specific 
signal can be established. 

II. WHY SHOULD GEO-NEUTRINO SPECTRA 
BE MEASURED? 

Geo-neutrinos are produced through pure [3 and 
processes: 

X -v X' + e + P e 
X -> X'*+e + 9 e 

\ 

X' + n + 7 . 

In order to determine the geo-neutrino decay spectra 
f{E v ) one has to know: 

(i) the feeding probabilities p n of the different energy 
states of the final nucleus; 

(ii) the shape of the neutrino spectrum for each tran- 
sition. 

Let us discuss in some detail the procedures and assump- 
tions used for deriving these quantities. 

Feeding probabilities are derived from measurements 
of the intensities E?> n of the gamma lines. These are 
corrected for internal conversion in order to derive the 
transition probabilities from level m to n: 

I m ' n = I™> n (l + a m ' n ) . (4) 

The internal conversion coefficients a m,n are obtained 
by theoretical calculations. In general they are of order 
10~ 2 , unless selection rules forbid or inhibit the gamma 
emission 2 . 

The feeding probabilities for the excited states are then 
obtained with a subtraction procedure as the difference 



1 A detector with N p free protons will collect a signal rate S = 
Np"^2i&iSi, where <I>; are the incoming fluxes of geo-neutrinos 
from the i-th beta decay in the chain and Sf are the corresponding 
specific signals. 

2 An important case in this respect is the EO transition of 214 Bi at 
1415.8 keV, which occurs essentially through internal conversion. 



between the intensities of outgoing and ingoing transi- 
tions: 

m<n m>n 

The feeding probability of the lowest state, po, is obtained 
with the same subtraction procedure: 

Po = l-Y, im, ° ■ ( 6 ) 

This procedure implies that all transitions to the ground 
state that are not observed, or taken into account, are 
included in the feeding probability to the lowest en- 
ergy state. In other words, po is indirectly determined, 
whereas it is of special interest for our purposes: (3 tran- 
sitions directly to the lowest energy state, the "pure /?" , 
produce the most energetic geo-neutrinos and thus give 
the largest contribution to the specific signal. 

For each transition, the shape of the neutrino spectrum 
is generally calculated assuming the well-known "univer- 
sal shape" distribution. This expression, see 0], corre- 
sponds to momentum independent nuclear matrix ele- 
ments (as for allowed transitions) and includes the effect 
of the bare Coulomb field of the nucleus through the rel- 
ativistic Fermi function. Electron screening and nuclear 
finite size effects are not considered. Note that this same 
"universal shape" expression is used even for the forbid- 
den transitions (see Tables U and [TTJ) , where momentum- 
dependent nuclear matrix elements can appear. 

These observations suggest that the feeding probabili- 
ties need to be confirmed by different experimental tech- 
niques and the electron decay spectrum need to be ex- 
perimentally tested. 



III. TOWARDS A DIRECT MEASUREMENT 
OF GEO-NEUTRINO DECAY SPECTRA 

When the nucleus X decays, whichever is the transi- 
tion involved, energy conservation provides a connection 
between the neutrino energy E v , the kinetic energy of the 
electron T e , and the total energy of the emitted gammas, 
Ej : 

Q = E v + T e + E 1 , (7) 

where Q = Mx — Mx 1 — M e is the Q-value for the decay. 
In order to measure the geo-neutrino spectrum, one needs 
a calorimetric detector which is capable of measuring the 
"visible" energy deposited together by electrons 3 and 
gammas, E v - ls = T e + E 1 . When measured decay events 
are displayed as a function of E V1S , by a mirror reflection 



3 Note that the energy deposited by conversion electrons is also 
included. 
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TABLE I: Effective transitions in the 238 U chain from In bold the reactions that give most of the signal. For each decay the 
table shows the probability, the maximal antineutrino energy, the intensity Ik, its error Alt, type and percentage contributions 
to the uranium geo- neutrino signal, and to the (U+ Th) geo- neutrino signal. For this last column it is assumed the chondritic 
ratio for the masses (Th/U = 3.9), which implies that 79% of the geo-neutrino signal comes from uranium. 
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TABLE II: Effective transitions in the Th chain from 9]. In bold the reaction that gives most of the signal. For each 
decay the table shows the probability, the maximal antineutrino energy, the intensity Ik, its error Alt, type and percentage 
contributions to the thorium geo-neutrino signal, and to the total (U+ Th) geo-neutrino signal. For this last column it is 
assumed the chondritic ratio for the masses (Th/U = 3.9), which implies that 21% of the geo-neutrino signal comes from 
thorium. 
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FIG. 1: The decay spectrum of 214 Bi as a function of the 
visible energy E v i s (full line) and of the neutrino energy E v 
(dashed line). 



one immediately obtains the number of events as a func- 
tion of neutrino energy, at E v = Q — E v - ls ; as an example, 
see figure Q] for the decay spectrum of 214 Bi. 

For such a measurement one needs a detector that can 
collect the energy lost by both electrons and gammas 
and that has a similar response to both particles. Es- 
sentially, this is a calorimetric measurement. In prin- 
ciple, it can be done with large bolometers which 
have very good energy resolution but long dead times. 
A sufficiently large liquid scintillator detector is suitable 
for such measurements. Although energy resolution is 
limited, nevertheless it can contain both electrons and 
gammas and a significant statistics can be collected in a 
reasonable time. 



There are some limitations that should be considered, 
when using a scintillator as a calorimeter. An ideal de- 
tector should provide the same response for gammas and 
electrons with equal energy, independently of the posi- 
tions where the particles are generated. In practice, how- 
ever: 

(i) even in a very large detector, the energy released 
as scintillation light from electrons and gammas of 
the same energy are not the same. This difference 
becomes marked at low energy, see Fig. [21 

(ii) Gammas can escape from a finite detector, thus 
releasing only a fraction of their energy. 

(iii) The number of photons collected by the detector 
can depend on the position where they have been 
produced (due to absorption, optical coverage, . . . ). 

All these effects can be taken into account by using 
calibration measurements, by selecting events that occur 
in the inner part of the detector (in order to minimize 
corrections due to escaping gammas) and with energy 
above a suitable threshold. The comparison between ex- 
perimental spectra and theoretical predictions has to be 
implemented by means of a Montecarlo simulation that 
accounts for the actual characteristics of the detector. 



IV. THE PROPOSED DETECTOR 

We propose to exploit the potential of the Counting 
Test Facility (CTF), which is operational and available 
in the underground I. N.F.N. Gran Sasso National Labo- 
ratory. 
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FIG. 2: Quenching factor for electrons and gammas, for a 
Birks coefficient ks = 1.95 ■ 10~ 3 cm MeV -1 , see appendix. 
The quenching factor has been normalized so that is one for 
I MeV electrons. 
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FIG. 3: Side view of the design of CTF. The vessel (labelled 
RI00 in this drawing) and shroud (R200) are shown, as well 
as the six rings of PMT's, the cylindrical tank, and the tubes 
used for filling and draining the vessel. The point PNTO is 
the nominal center of the sphere of PMT's and of the CTF 
vessel. Dimensions are given in cm. Courtesy of the Borexino 
collaboration. 



Like Borexino, the CTF design [13j is based on the 
principle of graded shielding, see figures [3] and 2J The ac- 
tive scintillation liquid in CTF is a four-ton mass of pseu- 
documene enclosed in a transparent nylon sphere, the 
CTF vessel. Outside this vessel there is a volume of ultra- 
pure water which is enclosed in a second nylon sphere, the 
so-called CTF radon shroud, intended to prevent radon 
transport with thermal fluxes from the outside zones of 
the detector. A set of inward-facing PMTs is arrayed 
outside the shroud. The entire apparatus, surrounded 
by another volume of water, is contained in a cylindri- 
cal stainless tank. The bottom surface of the tank holds 
16 upward- facing PMTs used to tag the muons passing 
through the detector by means of the their Cherenkov 
light in the water. 

The facility is equipped with a rod system, which can 
be used to insert a small, cylindrical quartz vial inside 
the CTF vessel. A suitable source, dissolved in the liq- 



FIG. 4: A picture of the CTF viewed from below. Courtesy 
of the Borexino collaboration. 

uid scintillator, can be placed in the vial. Electrons are 
stopped inside the vial and the scintillation light is prop- 
agated within CTF through the quartz (which is trans- 
parent to the near-UV wavelengths of scintillation light 
and has index of refraction close to that of the scintil- 
lator) whereas gamma conversion occurs inside the CTF 
inner vessel. The inward facing PMTs outside the shroud 
can thus detect light originating from both electrons and 
gammas. 

A Montecarlo code has been developed for CTF. It is 
described in [TJ] and in the Appendix with the adjust- 
ments for our specific task. 

V. WHAT HAS TO BE MEASURED? 

Geo-neutrinos with energy above the threshold for re- 
action |T]) arise only from the chains of 238 U and 232 Th. 

In particular, for 238 U only three nuclides ( 234 Pa, 
214 Bi, 210 T1) contribute to the geo-neutrino signal. The 
contribution from 210 T1 is negligible, due to its small oc- 
currence probability, and the uranium contribution to the 
geo-neutrino signal comes from five (3 decays: one from 
234 Pa and four from 214 Bi. TableUlists the effective tran- 
sitions, i.e. , those that can produce antineutrinos with 
energy above the threshold E$. In fact, 98% of the ura- 
nium signal arises from the two transitions to the ground 
state (in bold in Table Hi]) and an accuracy better than 
1% is achieved by adding the third one. 

232 Th decays into 208 Pb through a chain of six a- 
decays and four /3-decays. In secular equilibrium the 
complete network includes five /3-decaying nuclei. Only 
two nuclides ( 228 Ac and 212 Bi) yield antineutrinos with 
energy larger than 1.806 MeV. The thorium contribution 
to the geo-neutrino signal comes from three (3 decays: 
one from 212 Bi and two from 228 Ac (see Table HI)) . In 
fact, 94% of the thorium signal arises from the transition 
to the ground state of 212 Po (in bold in Table HT|) . 

We remind that, assuming the chondritic ratio for 
the global uranium and thorium mass abundances, 
a(Th)/a(U) = 3.9, one expects that geo-neutrinos from 
uranium (thorium) contribute about 80% (20%) of the 
total U + Th geo-neutrino signal. 

In summary: 
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(a) 98% of uranium geo-neutrino signal comes from 
just two transitions, one from 214 Bi and the other 
from 234 Pa. They provide 77% of the expected to- 
tal U + Th signal. 

(b) A single decay of 212 Bi accounts for 94% of the 
thorium signal. It provides 20% of the expected 
U + Th signal. 

Just three transitions have to be investigated experimen- 
tally. In this respect, the following considerations can be 
useful: 

(a) 222 Rn (ti/2 = 3.824 days) can be easily dissolved in 
the scintillator and the decay of 214 Bi is uniquely 
identified by the subsequent decay of 214 Po (r X j 2 = 
164.3^s). 

(b) By dissolving 238 U in the scintillator, one can de- 
tect the beta decay of 234 Pa (superimposed, how- 
ever, with that of 234 Th). The subsequent decays 
of the chain are effectively blocked by the long half- 
live of 234 U (t 1/2 = 2.455 ■ 10 5 yr). 

(c) For the investigation of 212 Bi decay one has to start 
with a 224 Ra source (ti/ 2 — 3.66 days) or with 
a 232 Th source. The decay of 212 Bi can be eas- 
ily identified by the subsequent a decay of 212 Po 
(t 1/2 = 299 ns). 




FIG. 5: Data and best fit. Events have been grouped into 
65 bins from to 3.4 MeV. Event numbers together with 
statistical error (vertical bar) and bin size (horizontal bar) are 
presented as a function of bin number, see text. We only fit 
data in bins from 3 to 65. The continuous curve shows the best 
fit when three parameters (po, light yield, and normalization) 
are left as free. 
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FIG. 6: A simplified decay scheme of 214 Bi and the contri- 
bution of the various levels to the total geo-neutrino signal. 



VI. RESULTS FROM A DIFFUSE RN SOURCE 

In order to test the method that we are proposing, 
we have used data from a sizeable, though limited, radon 
contamination of CTF, which occurred in the early phase 
of operation of the detector. In the full volume of CTF, 
we selected the candidates /3-decay from 214 Bi by the 
distinctive subsequent 21 Po a decay, which occurs with 
a mean time-delay of 237 (is (the so-called BiPo events). 

The selection of the analyzed events, see the data 
points in Fig. [5l is described in the following two sub- 
sections IVFAl and IVTBl 

We recall that most of the contribution to the geo- 
neutrino signal arises from the transition to the lowest 
energy state (0) of 214 Po, see Fig. [SI Our analysis aims: 

(1) to determine the probability po of populating the 
lowest energy state (assuming the universal "al- 
lowed" shape) from the CTF data. 

(2) To determine whether the spectrum of the pure 
beta transition (that to the lowest state) is de- 
formed with respect to the universal allowed shape. 

(3) To discuss the implications of this study on the spe- 
cific geo-neutrino signal s( 214 Bi), given by Eq. ^j. 



A. Data selection and backgrounds 

The main selection criterion is that the coincidence 
time between consecutive signals, provided by the 
prompt (5 decay of 214 Bi and the delayed a decay of 
214 Po, be 2/is < Ai 1/2 < 602^s. The lower limit elimi- 
nates fast coincidences from the corresponding decays in 
the 232 Th chain, whereas the higher limit has been chosen 
to keep random coincidences under 1%, while preserving 
high statistics. The selected events are 4.54 • 10 5 . 

Then we require that the energy deposited by the first 
signal is E\ < 3.9 MeV, taking into account the Q- 
value (Q = 3.27 MeV) and the energy resolution, about 
0.2 MeV at these energies. This cut removes random co- 
incidences, while keeping the Bi-Po events (acceptance 
almost 100%). 

The reconstructed radial positions 4 of the two signals 
n and r 2 , are used to impose the three conditions: n < 
2 m, r 2 < 2 m, and \r\ — r 2 | < 2.5 m. These very weak 
cuts have total acceptance efficiency about 100%, while 
removing random coincidences. 

We then impose that the energy of the second sig- 



4 A typical spatial lcr resolution is of order 10 cm at 1 MeV. 
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nal, on the electron energy scale, is 0.56 MeV < E2 < 
1.1 MeV (note that the 7.9 MeV alpha particles from 
214 Po decay are quenched by a factor about 11), in order 
to reduce low-energy a's of the Rn chain and random 
coincidences. The acceptance efficiency is 98.7%. 

CTF has goo d pulse-shape discrimination between a 
and (3 events [15j . To avoid contamination by high-energy 
a's, we add cuts on a suitable a/ (3 discrimination param- 
eter. The combined acceptance efficiency is 99.4% and we 
are left with 4.46 • 10 5 events. 

The remaining background is estimated by applying 
the same sets of cuts, with the coincidence time window 
2000/iS < Ai x / 2 < 8000/is. In this way, we estimate that 
contamination of random coincidences is about 0.7%. 



B. Data analyses and fiducial region 

Taking into account the number of events and the esti- 
mated energy resolution AE ss 80 keV ■ (E/MeV) 1 ^ 2 , we 
grouped events in bins of about 50-keV 5 and analyzed 
the 63 bins from 0.1 MeV up to 3.4 MeV. 

In order to reduce systematic effects due to 7's that 
are only partially contained and due to deviations from 
spherical symmetry, one should select events near the 
detector's center. On the other hand, statistics improves 
with increasing volume. We found that a good compro- 
mise is to consider events such that the a's reconstructed 
positions are within a sphere of 42 cm around the CTF 
center. This sphere, after the cuts discussed above, con- 
tains 3.14 • 10 4 candidate decays. 

The theoretical spectra have been produced with the 
CTF code, described in [lj] and with the specific adjust- 
ments presented in the Appendix. 



C. Feeding probability of the lowest state 

First we shall assume that the neutrino energy distri- 
bution f{E v ) is given as a sum of "universal" functions, 
i.e. , 

PnF univ 

(E V ,Q-E n ) , (8) 

n 

where E n is the energy of the n-th level {Eq = 0), i.e. , 
the maximal energy which can be taken by the neutrino 
is Q — E n , and the functions F un - lv {E v ,Q — E n ) are each 
normalized to unity. The electron kinetic- energy distri- 
bution is: 

<j>(T e ) = 5> n $ im iv(T;, Q - En) . (9) 
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FIG. 7: Data and best fit: residuals. The values shown are 
the residual relative to figure [5] data - best-fit value divided 
by the square root of data. The best fit is obtained with three 
parameters (po, light yield, and normalization) left as free for 
data in bins from 3 to 65. 

The universal distributions for neutrinos and electrons 
are related by: 

$univ (T e ,Q — E n ) = F univ [E u = Q — E n — T e ,Q — E n ) . 

(10) 

The populations of the 82 excited 214 Po states are fixed at 
the values given in the Table of Isotopes (Tol) [l6| , apart 
for a common normalization factor, such that the total 
population of these states is (1— po). This assumption for 
the excited states means that the relative intensities of 
the 7 transition lines in the decay are exactly determined. 

We fitted the data with Montecarlo generated spectra 
leaving as free parameters: 

(1) po, the feeding probability of the lowest state; 

(2) the Light Yield L, defined as the number of pho- 
toelectrons that would be collected by 100 photo- 
multipliers for an electron depositing 1 MeV at the 
center of CTF; 

(3) the normalization, i.e. the number of reconstructed 
candidates (which should be equal to the number 
of candidates). 

The best fit function and the residuals are shown in 
Figs. E] and [3 At the minimum x 2 /d.o.f. = 61.7/(63-3), 
the Light Yield L = 321 p.e./MeV 6 , and the normaliza- 
tion factor is 0.998. The best fit value is po = 0.177 with 
a statistical la error of ±0.004. 

Systematic uncertainties arise from limitations of the 
Montecarlo simulation with respect to the real detector 
(see Appendix). We found that the largest uncertainties 
originate from the imperfect spherical symmetry of the 
detector arising because of the deformations of the inner 



5 The energy scale depends on the best-fit value of the light yield. 
At the best fist the bin size is 52 keV. 



Note that this Light Yield is within 3% the one determined at 
much lower energy by a fit to 14 C events. 
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vessel (IV) and to non-spherical distribution of the ac- 
tive PMTs. The nature of these systematic errors makes 
them more important for events at large radii or involv- 
ing high-energy gammas, which can deposit energy far 
away from the point of origin. Therefore, we estimated 
the effect of these errors on our measurements by ana- 
lyzing subgroup of events with different distance r from 
the center of the detector. Results show a consistent be- 
havior as function of r. Effects of the uncertainties on 
the quenching parameters, subtraction of random coin- 
cidences, selection of the energy window, and choice of 
the size of the energy binning have also been considered. 
In conclusion the total systematic error is estimated as 



h0.003 t - , 
-o.ooi: ho Ln ai. 



(CTF) p = 0.177 ± 0.004(stat)^;^(sys) 



(11) 



This value is consistent with that reported in Tol [la ]. 
Po(ToI) = 0.182 ±0.006. 

D. Shape factor for the pure beta transition 

Next we release the assumption that the spectrum for 
the transition to the ground state has the "universal 
shape". The electron energy distribution is assumed to 
be: 

4>{T e ) =p $(T e ) +^p„$univ(T e ,Q- E n ) (12) 



n>0 



where 



$(T e ) = $ u „iv(T e ,Q) 1 + 2/ 



Te ~ (TV 



(Te) 



(13) 



and the average energy (T e ) is calculated over 
$u„iv(T e ,Q). 

The dimcnsionless "shape parameter" y describes thus 
a deviation from the universal formula. Note that this 
simple parameterization does not change the normaliza- 
tion of the distribution; it only changes its shape. Other 
parameterizations are of course possible: at the level of 
the accuracy of this preliminary study, they would not 
change our conclusions. 

Present data do not allow to independently determine 
Po, pi, and y. We consider therefore as inputs the values 
given in Tol, p = 0.182 ± 0.006 and p x = 0-017 ± 0.006, 
with errors assumed uncorrelated, and leave only y as 
an unconstrained parameter. The resulting % 2 is shown 
in Fig. H as a function of y. At the best fit we find 
X 2 /d.o.f. = 51.6/(65 - 5), p = 0.177, Pl = 0.008 and: 



y : 



-0.11 ±0.06(stat) 



(14) 



The universal spectrum, y = 0, has a x 2 larger by 5.9 
with respect to the minimum: the statistical evidence for 
a deformed spectrum is about 2.4<t. If both p and p\ are 
left completely unconstrained, one again finds that the 
best fit value for y is —0.11, but with a larger la interval: 




FIG. 8: The \ 2 as a function of the deformation parameter 
y. The number of degrees of freedom is (65 — 5). 



—0.53 < y < —0.09. However, spectral deformation is 
very sensitive to the lowest part of the visible energy 
and, therefore, we expect larger systematic uncertainties 
than the one estimated in the case of po- In fact, the 
analysis of subgroup of events with different distances 
from the center of the detector, as we did for the feeding 
probabilities, results in an estimated systematic error for 
y of the same size of its deviation from zero, i.e. ss 0.10. 

Our present result only shows that the method is sensi- 
tive to the form of the spectrum and has the potentiality 
of detecting spectral deformations. However, interesting 
results can be obtained by achievable reductions of sta- 
tistical and systematic errors. 



E. Implications for the specific signal 

The geo-neutrino signal, s( Bi), can be written as 
the sum of two contributions: 



s( 214 Bi) =p (a) +pi(a)i 



(15) 



where the cross section of reaction (fTJ) is averaged over 
the neutrino energy distribution. Assuming universal 
shape, i.e. , (<r) = / dE l ,a(E l/ )F nn i v (E v , Q - E n ), the 
cross sections are (a)o = 7.76 and (a)i = 2.825 in 
units of 10 -44 cm 2 , with errors of order half of percent. 
Clearly, the largest contribution to the geo-neutrino sig- 
nal is given by the first term in Eq. (fTS"]) . so that the 
relative error is practically the one on pq. From the pre- 
vious analysis we find: 

(CTF) (16) 
s( 214 Bi) = (l.42±0.03(stat)±g;ggi(sys)) • 10~ 44 cm 2 . 

This should be compared with the result derived using 
the Table of Isotopes: 



(Tol) 

s( 214 Bi) = (1.46 ± 0.05(stat)) ■ 10~ 44 cm 



(17) 
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If spectral distortion is allowed in the form of Eq. (fT3]l . 
then (cr)o becomes 

(<j)o -> (a)o+J/(Aa) , (18) 

where (ct)q = 7.76 and (Acr}o = -4.52. 

If we substitute the value of y in Eq. (|14p and the 
corresponding values for p = 0.177 and pi = 0.008, we 
find: 

s( 214 Bi) = (1.48 ± O.Ol(stat) ± 0.03(sys)) • 10~ 44 cm 2 . 

(19) 

Note that, if we leave completely unconstrained the 
shape, po, pi, and y, the effects of changes of shape and of 
p on the signal are anti-correlated: if the spectrum is de- 
formed so that there are more (less) low-energy electrons, 
the corresponding best fit value for po is lower (higher). 
For instance using the present parameterization, we can 
let y span from -0.64 to 0.13, finding the corresponding 
best fit values for po e p\ with no constraint (we disre- 
gard the fact that these best fit values have often too 
large x 2 ) : while the values of po go from 0.13 to 0.20, the 
signal changes only by about ±2%. In other words, the 
resulting signal is weakly dependent on the shape factor. 

VII. CONCLUDING REMARKS 

So far, we have estimated the 214 Bi geo-neutrino spe- 
cific signal by using CTF data resulting from a limited 
radon contamination. Our estimate has a comparable er- 
ror with that derived from Tol, see Eqs. (fi"6|) and (TiT]) . 
We remark, however, that our method has two advan- 
tages: 

(i) the pure beta transition can be detected in CTF 
and its probability can be measured directly 
(whereas in the study of gamma lines alone its ex- 
istence was inferred from the fact that the gamma 
counts did not match with the expected number of 
decays and the probability was evaluated from this 
mismatch) ; 

(ii) one can check the validity of the universal shape 
approximation for the most important decay mode. 

A dedicated experiment makes sense, with the aim of 
reducing the relative statistical error Apo /po to the level 
of the relative error on the cross section, A(cr)o/(cr)o ~ 
0.5%. This requires a statistics larger by a factor of about 
20, or some 6 • 10 5 selected events. At the same time, one 
has to reduce systematic errors on the correspondence 
between measured light and released energy: the largest 
improvement should be obtained by concentrating the 
source near the center of the detector. 

Our preliminary results are of encouragement towards 
a series of dedicated measurements, with suitable sources 
dissolved in the liquid scintillator and placed in a vial 
near the CTF center. As an example, a 222 Rn source 
with an initial activity of 5 Bq would be tolerable for CTF 



and it would produce some 2 • 10 6 decays in 11 days (2 
lifetimes). Electrons are stopped inside the vial and the 
scintillation light is propagated within CTF through the 
quartz, whereas gamma conversion occurs inside the CTF 
inner vessel. The inward facing PMTs outside the shroud 
can thus detect light originating from both electrons and 
gammas. Along these lines, one can get a better estimate 
of the specific signal of 214 Bi, provide measurements of 
the other signals, s( 234 Pa) and s( 212 Bi), relevant for geo- 
neutrino studies, and, more generally, measure the energy 
spectra of neutrinos from long lived heavy nuclei. 
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APPENDIX: MONTECARLO SIMULATION OF 
CTF DETECTOR 

Because light propagation in a large-volume scintilla- 
tor detector involves complex mechanisms, the precise 
modeling of the detector response requires that many 
phenomena are taken into account. Among the most rel- 
evant issues, it is worth mentioning the wavelength de- 
pendence of the processes involved in light propagation, 
the reflection/refraction at the scintillator/ water inter- 
face, and the light reflection on the concentrators. The 
need to follow each of some 12000 photons emitted per 
1-MeV electron event makes tracing MC code very slow. 

A fast and reliable code has been developed for the 
CTF detector and it is briefly described in this appendix; 
more details can be found in [l4j. The code takes ad- 
vantage of using average parameters (such as light yield, 
energy resolution, and spatial reconstruction precision) 
obtained analyzing the detector's data. Optimal sets 
of data for calibrating, tuning and testing the code are 
the /?-decay spectrum from 14 C and the easily identifi- 
able a's from the radon-chain decays. The code has two 
parts: the electron-gamma shower simulation (EG code) 
and the simulation of the registered charge and position 
(REG code). 

The EG code generates events at a random position 
with random initial direction (for 7's) and follows the 
gamma-electron shower using the EGS-4 code [l7| . The 
low-energy electrons and alphas are not propagated in 
the program and are considered to be point-like sources, 
located at the initial coordinates. The mean registered 
charge corresponding to the electron kinetic energy T e is 
calculated with 

Q e (r)=A-T e -f(k B ,T e )f R (r) , (A.l) 
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where fa(r) is a radial factor, which takes into ac- 
count the dependence of the registered charge on the dis- 
tance from the detector's center, and f(kB,T e ) is the 
ionization quenching factor for electrons; the normal- 
ization of these two factors has been chosen such that 
f R (0) = f(k B ,T e = 1 MeV) = 1. The method used to 
obtain fn{r) consists in studying the response for mono- 
energetic alpha's as a function of their radial position and 
it is described in (l8| . 

For the PC-based scintillator (PC+PPO 1.5 g/1), the 
quenching factor k B = (1-7 ± 0.1) • 10 3 cm MeV -1 was 
found to satisfy experimental data . This value agrees 
with the fit to high-statistics (3 spectrum of 14 C. The a 
particles from 214 Po decay, which tag the Bi-Po events 
considered in our study, have energy of 7.69 MeV and are 
quenched to an equivalent /3-energy (produce the same 
amount of light in the scintillator of an electron) of 751 ± 
7 keV. In the set of data selected for the present work 
quenching was higher (and the light yield lower) than in 
Ref. due to the presence of oxygen in the scintillator 
(the radon originated from atmospheric air). In fact, 
the 7.69 MeV alphas are found at the lower equivalent 
f} energy, E = 643 keV, for events selected around the 
detector's center. The ratio of two energies can be used 
to scale the quenching k B factor: the adopted value is 
k B = 0.0195 cm MeV" 1 . 

The 7's are propagated using the EGS-4 code. As 
soon as the i-th electron of energy T ei appears inside the 
scintillator, the corresponding fraction of total registered 
charge is calculated: 

AQi = A-T ei -f(kB,T ei )f R (ri) . (A.2) 

The total mean collected charge is defined, when the 7 is 
discarded by the EG code, as the sum of the individual 
deposits: 

Q 7 = ^AQ 4 . (A.3) 

i 

The weighted position is assigned to the final 7: 

_ E» A Q» • r » (A4) 



where AQi is the charge deposited by the i-th electron 
at the position r,. 



Once the position and deposited charge of the event 
have been generated by the EG code as described above, 
the second part of the code (REG) randomly generates 
the corresponding number of photoelectrons registered at 
each PMT; it takes into account the proper geometrical 
factor and assumes Poissonian distributions of photoelec- 
trons number at each PMT. 



Finally, the energy-dependent radial reconstruction is 
simulated. The reconstruction precision is assumed to 
be defined by the number of PMTs fired and dependent 
only on the distance from the detector's center (spherical 
symmetry) . These two assumptions have been confirmed 
by measurements using artificial radon sources inserted 
in the CTF detector [lEHS]- 



The main source of systematic errors in our study is the 
departure from spherical symmetry of the detector due 
to deformations of the inner vessel (IV) and non-uniform 
distribution of active PMT on the spherical surface sur- 
rounding the scintillator. The IV, a 500/^m thick nylon 
bag containing four tons of low density (« 0.88 g/cm 3 ) 
scintillator, is immersed in water. The buoyancy forces 
are compensated by supporting strings, but the deforma- 
tions are not measured precisely and are not accounted 
for in the MC modeling. The maximum radial deviations 
from the ideal sphere can be as big as 5 -j- 10 cm, though 
in average the radius of the sphere is R = 100 cm. An- 
other source of systematics of the same nature (absent 
in an ideal spherical detector) is the position dependence 
of the light-collection efficiency function, fn(r), which is 
assumed to depend only on the distance from the center 
and not on all three coordinates. The nature of these sys- 
tematic errors makes them more important for events at 
large radii or involving high-energy gammas, which can 
deposit energy far away from the point of origin. 
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